Monte Carlo Study of an Extended 3-State Potts Model on the Triangular Lattice 
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By introducing a chiral term into the Hamiltonian of the 3-state Potts model on a triangular lattice 
additional symmetries are achieved between the clockwise and anticlockwise antiferromagnetic states 
and the ferromagnetic state. This model is investigated using Monte Carlo methods. We investigate 
the full phase diagram and find evidence for a line of tricritical points separating the ferromagnetic 
and antiferromagnetic phases. 
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PACS numbers: 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

The-, 3-state nearest neighbour ferromagnetic Potts 
modetl on the triangular lattice has a second order phase 
transition at a finite temperature which can be deter- 
mined from duality and a star-triangle relation, together 
with the assumption of a single critical pointa. The criti- 
cal exponents are known exactly to be a = | , (3 ■ 



hi 



^ and v = |. However, the anti- ferromagnetic model 
displays a completely different behaviour. It is believed 
to have a weak first ,-osder transition to a six- fold de- 
generate ground state&oEI In this rftaper we discuss an 
extension of the 3-state Potts modelErQ to include a chiral 
term which clearly exposes the first order nature of this 
latter transition. 



II. THE MODEL 



The 3-state Potts model is usually defined as 



(1) 



i<j 



where the sum is over nearest neighbours and the vari- 
ables <ii can have three values A, B, C at each site. For 
the triangular lattice, it is convenient to divide the lattice 
into three equivalent interpenetrating sublattices corre- 
sponding to the three sites on each triangle. The Hamil- 
tonian can then be written as a sum over all triangles 
A if nearest neighbour bonds are shared between neigh- 
bouring triangles: 



H = J2Ha 



J 



where <j\ , 02 and 173 correspond to the three sublattices 
on each triangle. 

We now extend this model by adding two additional 
three-spin terms to each triangle 



Ha = -r + (r-^)(^ lff2 
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(3) 



where e 



is the Levi-Civita symbol. The term in- 



volving £1 is a chiral term which distinguishes between 
cyclic and anti-cyclic permutations of the three states 
on a triangle. The terms involving T are equivalent to 
adding the term — re^ icr , 3(T3 to each triangle. Each trian- 
gle has 27 states corresponding to four different energies 
as indicated in Table ||. 

This extended Potts model was first introduced and 
studied using real space renormalization group methods 
by Young and LavisO. Theirj'esults agreed with previous 
work by Schick and GriffithsEl on the conventional 3-state 
Potts model which predicted a second order phase tran- 
sition for the antiferromaflnet on this lattice. However, 
more recent Monte Carlcfffl results predict a weak first 
order transition as well ft&-,a tricritical point located in 
the ferromagnetic regiorJj'El A related self-dual model 
which only has the conventional three spin terip-on the 
upward triangles also exhibits a tricritical point)13. 

In the the present work, we consider the extended Potts 
space which is described by the three independent cou- 
pling parameters J, T and f2. The relationship of our pa- 
rameters J, r and fl tOjlhe parameters K and M defined 
by Schick and Griffithaj ancLthe additional parameter P 
defined by Young and LavisEl is 



(2) 



K = (J-2T)/T 
M = (3J-2T)/T 
P = 2CL/T 



(4) 



2 




where T is the temperature. 



III. SYMMETRIES 

The model has a symmetry under the following trans- 
formations: If we label the three states as A, B, C and the 
three sublattices as 1, 2, 3 as in Table |, then the transfor- 
mation which leaves the states on sublattice 1 unchanged 
but performs a cyclic permutation of the states on sub- 
lattice 2 and an anti-cyclic permutation of the states on 
sublattice 3 maps the configurations on each triangle as 
follows: 

ei — ► e 3 , e 2 — > e 2 , e 3 — > e 4 , e 4 — > ex (5) 

and the partition function is invariant under all permu- 
tations of the parameters ex, e 3 and e 4 . 

A convenient set of parameter^ to label special points 
and lines in parameter space area 

e 3 - ex 

a = 

e2 - ex 



The parameters J, and T are related to a and 6 as 

17 b — a 

1 = 2 

r 3 — a — b 

I = —5— (7) 

Figure 1 shows the possible ground states and the lo- 
cation of models corresponding to special values of (a, b) . 
Under the symmetries described above, all cases can be 
mapped into the upper triangular region. 

The point F(3/2,3/2) is the usual 3-state nearest 
neighbour ferromagnetic Potts model which corresponds 
to il — 0, r = and AF is the corresponding anti- 
ferromagnetic model which maps under the symmetries 
of the extended model to the point (0,3). The point 
T(l, 1) corresponds to J = 2r, fi = and is equivalent to 
a 3-state Potts model with only the usual three spin term 
non-zero and coupling J. The point Q (1, 2) is an image 
of the model with only the chiral term f2 non-zero. All 
ordered phases meet at the point G (0,0). 
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FIG. 1: The possible ground states and the location of models 
corresponding to special values of (a, b). Under the symme- 
tries described above, all cases can be mapped into the upper 
triangular region. 

IV. RESULTS 

The extended Potts model studied here, can also be 
represented in terms of a spin-1 Ising modelQ. Identifying 
the 3 states A, B, C as the spin states Si — ±1, at each 
site i, we choose the following order parameter 

m = |mx| + |m2| + |m 3 | (8) 

where m a , a = 1, 2, 3 are the sublattice spin magnetiza- 
tions. At high temperatures, all three states have equal 
weights and the sublattice magnetizations are zero apart 
from finite size effects. At low temperatures, the sublat- 
tices become ordered into one of the three possible states. 
In the ferromagnetic phase, the sublattices have the same 
state. However in the anitferromagnetic regions, there 
can be either a cyclic or anti-cyclic arrangement. We al- 
low for the permutation symmetry of each of these phases 
by summing m over the permutations. In this way m has 
the value unity in each of the ordered phases at zero tem- 
perature. We also calculate the susceptibility x an d the 
Binder cumulant Um associated with fluctuations in this 
order parameter. 

We have used a standard Metropolis Monte Carlo al- 
gorithm to study lattices with N — L 2 sites for values 
of L ranging from 18 to 96 along the line T = joining 
(3/2,3/2) to (0,3), followed by the line Q + T = 3J/2 
from (0,3) to (0,0) and finally back along the line £1 = 
from (0,0) to (3/2,3/2). We present our results for several 
points along these lines: 

Point Q(l,2) 

This point is equivalent under the permutation symme- 
try to the case where only the parameter fl is non-zero. 
The system is expected to exhibit a second order phase 
transitionO belonging to the same universality class as 
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FIG. 2: Binder cumulant as a function of T/J for sizes L — 
18,24,30,36,42,48 in the case of a = 1,6 = 2. The curves 
cross at T c /J = 1.5623 indicating a continuous transition. 



the ferromagnetic 3-state Potts model. Our Monte Carlo 
results are consistent with this behaviour. We have cal- 
culated the order parameter (to), the susceptibility (x) 
and the Binder cumulant {Urn) for lattice sizes ranging 
from L = 12 to L = 48 where N = L 2 is the num- 
ber of sites. Figure 2 shows the Binder cumulant Um 
from which the crossing point of the curves for different 
L yields an estimate of the critical temperature. Finite 
size scaling is then used to collapse the data for both the 
order parameter and the susceptibility (figure 3) in order 
to determine the critical exponents j3, 7 and v. We find 
13 = 0.11±0.01,7 = 1.44±0.01 and v = 0.82±0.06 which 
are consistent with this model belonging to the same uni- 
versality class as the ferromagnetic 3-state Potts model. 

Point AF (0,3) 

This point is equivalent under the permutation symme- 
try to the antiferromagnetic 3-state Potts model and is 
believed to exhibit a first order thermal phase transition. 
Our Monte Carlo results agree with those of Alder et. 
al.tJ. Figure 4a shows a histogram for the energy per site, 



e = E/N, at T c in the case of an 84 x 84 lattice and figure 
4b shows how the energies of the two peaks depend on 
lattice size. The results indicate a finite jump in the en- 
ergy in the limit of an infinite lattice. The temperature 
and energy scales at the point (0, 3) are twice as large as 
at the equivalent AF point in the lower part of figure 1. 

We have also studied this latter point directly with 
O = r = and J < 0. A histogram of the chirality order 
parameter e ai <j 2 a- 3 for an 18 x 18 lattice clearly indicates 
the first order nature of the transition. Figure 5 shows 
the histogram at a few temperatures near T c which, as 
mentioned above, has one half the value of the equiva- 
lent point at (0,3). At temperatures above T c , there is 
a single peak in the distribution at e — but, as the 
temperature is lowered, two additional side peaks appear 




FIG. 3: Finite size scaling plots for the order parameter and 
susceptibility for different lattice sizes L and reduced tem- 
peratures t — T ~ Ta in the case a = 1,6 = 2. a) The order 

J c 

parameter m near T c yields the values p — 0.11 ± 0.01 and 
v = 0.82 ± 0.06. b) The susceptibility \ yields the values 
7= 1.44 ±0.01 and v = 0.82 ±0.06. 



at finite non-zero values of e. At T c all three peaks have 
equal weight and below T c the two side peaks dominate. 
This behaviour clearly indicates that the phase transition 
is first order. A continuous transition would correspond 
to the central peak splitting into two peaks which move 
continuously away from zero. 

Point (0,1) 

This point lies on the first order surface separating the 
AAA and ABC ground states. We have calculated his- 
tograms for both the energy and the order parameter m 
for various lattice sizes L at the corresponding T c . Both 
distributions show a double peaked structure as T c is ap- 
proached. Figure 6 shows how the peak energies and 
magnetizations vary with L. Again we find a first order 
thermal phase transition which is stronger than at the 
AF point. 
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FIG. 4: a) Histogram for the energy per site, e = E/N, at 
T c /J = 1.256 in the case a = 0,fe = 3 which is equivalent 
to the 3-state Potts antiferromagnet. b) Dependence of the 
histogram peak energies on lattice size. Extrapolation to the 
large L limit indicates a finite jump in the energy at T c . 
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FIG. 5: Histogram for the chirality e at various temperatures 
near T c at the AF point with J < 0. 
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FIG. 6: Dependence of the histogram peak locations on 
lattice size for the case a = 0, b = 1. a) Dependence of the 
histogram peak energies on lattice size, b) Dependence of the 
histogram peak magnetizations on lattice size. Extrapolation 
to the large L limit indicates a finite jump in the both the 
energy and the order parameter at T c . 



Point G(0,0) 

This point corresponds to a meeting of the three first 
order planes separating the AAA, ABC and ACB ground 
states. Figure 7 show a magnetization histogram for a 
48 x 48 lattice at T c . The jump in the magnetization is 
larger than at the point (0,1)- We have also studied the 
behaviour at a large value of b = 10 along the line a = 
and again find a first order transition which is weaker 
than at the AF point. 

We have also investigated the possibility of a tricritical 
point along the line a = b using the histograms. Figure 
8 shows the energy discontinuity as a function of a. The 
results were obtained by calculating histograms for var- 
ious lattice sizes and extrapolating to the large L limit. 
The energy jump decreases as we move upward from the 
point G (0, 0) . A tricritical point exists very close to 
a = b = 0.35 which is in between the values reported by 
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FIG. 7: Magnetization histogram at the point G (0,0) for a 
48 x 48 lattice at T c . The peaks of the distribution are more 
widely separated than at the point (0,1). 
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FIG. 8: The energy jump as a function of a along the line 
a — b. A tricritical point exists near a = b = 0.35. 



the seriescl and previous Monte Carlo resultsa. 

Figure 9 shows the phase diagram along the closed path 
(3/2,3/2) -> (0,3) -> (0,0) -> (3/2,3/2) obtained using 
our Monte Carlo method. Along the first section of the 
path, the second order transition changes to a first order 
transition at the point labelled T which is quite close to 
the AF point. As we then move towards the point G (0, 0) 
along the line a = 0, the transition remains first order 
with the transition temperature decreasing and the size 
of the energy and order parameter discontinuities increas- 
ing. Finally, along the section of the path moving back 
to the point F the first order transition again changes to 
a second order transition near the point a = b = .35. The 
points labelled by T are tricritical points. 

Figure 10 shows a schematic phase diagram in the a — b 
plane where shaded areas correspond to first order tran- 
sitions. The shaded region narrows appreciably near the 



FIG. 9: The phase diagram along the closed path V — — > 
Q + T = 3J/2 — > = obtained using a heat-bath Monte 
Carlo method. Solid lines indicate second order transitions 
and dashed lines indicate first order transitions. Tricritical 
points are indicated by the letter T. 




FIG. 10: A schematic picture of the phase diagram in the ex- 
tended Potts space. The shaded area is a region of first order 
thermal phase transitions. The region narrows substantially 
near the AF point indicating a weak first order transition. 



AF point which is consistent with a weak first order tran- 
sition. 



V. SUMMARY 

We have studied an extended 3-state Potts model 
on the triangular lattice using Monte Carlo methods. 
The extended model includes a chiral term which distin- 
guishes between cyclic and anti-cyclic configurations on 
each triangle and clearly exposes the first order nature of 
the nearest-neighbour anti-ferromagnet. 

In the regions corresponding to continuous transitions, 
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the model belongs to the 3-state ferromagnetic Potts 
model universality class. We have not been able to locate 
the tricritical line with sufficient precision to determine 
the tricritical exponents. 
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